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ABSTRACT 

The  direction  of  sand  movement  along  the  Carmel  River  State  Beach 
was  qualitatively  determined  by  diving  observations,  a  bathymetric  survey, 
wave  refraction  diagrams  and  a  sediment  size  analysis  of  18  samples. 

The  primary  source  of  sediments  for  the  beach  appears  to  be  the 
Carmel  River  which  flows  only  seasonally.  Sedimentary  material  is  intro- 
duced into  the  bay  after  winter  precipitation  provides  a  sufficient 
amount  of  run-off  to  warrant  the  opening  of  the  river  mouth  by  bulldozer. 

The  fine  sedimentary  material  is  lost  offshore  and  the  coarser  ma- 
terial is  either  redeposited  on  the  beach  or  is  carried  south  with  the 
littoral  drift  and  deposited  at  the  nodal  point  in  the  sand  transport 
pattern.   This  node  is  located  on  the  northern  edge  of  the  head  of  the 
Carmel  Submarine  Canyon . 

Winter  storms  probably  induce  slumping  or  gravity  sliding  and  much 
of  the  material  is  carried  to  deeper  water  by  the  canyon. 
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I.      INTRODUCTION 

A.  OBJECTIVE 

This  study  was  undertaken  in  order  to  determine  the  direction  of 
sand  movement  along  the  Carmel  River  State  Beach,  Carmel,  California, 
and  provide  an  insight  into  the  methods  by  which  this  sand  is  introduced 
into  the  Carmel  Submarine  Canyon.  In  order  to  accomplish  these  objec- 
tives, the  area  of  study  was  examined  by  the  author  using  SCUBA,  a  bathy- 
metric  survey  was  conducted,  wave  refraction  diagrams  were  constructed, 
an  attempt  was  made  to  measure  bottom  currents,  and  sediment  samples  were 
taken  for  textural  analysis. 

B.  AREA  DESCRIPTION 

Carmel  Bay,  located  about  9  km  south  of  Monterey  Bay  on  the  central 
California  coast  (Fig.  1 ),  is  an  embayment  that  contains  the  head  of 
the  Carmel  Submarine  Canyon  (Fig.  2).   It  is  a  very  open,  rectangular- 
shaped  bay  with  dimensions  of  k-3   km  by  3.6   km  in  the  north-south  direc- 
tion and  east-west  direction,  respectively. 

Neighboring  communities  include  Pebble  Beach  and  Carmel,  California. 
At  the  southern  end  of  the  bay  is  Point  Lobos  State  Reserve.  Two  water- 
courses, the  Carmel  River  and  San  Jose  Creek,  drain  a  land  area  of 
approximately  670  km  and  terminate  near  the  head  of  the  Carmel  Canyon. 
[California  State  Department  of  Water  Resources,  1 969 ] - 

The  study  area  (Fig.  3  and  h)   is  bounded  by  Carmel  (Abalone)  Point 
to  the  north  and  San  Jose  Creek  to  the  south.  A  depth  of  75  ft  was 
arbitrarily  chosen  for  the  western  boundary  and  the  beach  face  formed 
the  eastern  boundary. 
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Figure   1 .      Location  of  Carmel  Bay  Showing  Watercourses  that  Drain 
into  It . 


Figure  2.  Carmel  Bay  Bathymetry  (Contour  Interval  -  10  fm) 
(After  Zardeskas,  1971). 
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Figure  3»     Area  of  Study. 


11 


Figure  k-      Aerial  Photograph  of  Southern  End  of 
Carmel  Bay  Showing  the  Area  of  Study. 


Figure  $.  Bulldozer  Opening  the  Mouth 
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The  primary  source  of  sediments  for  the  study  area  was  assumed  to  be 
the  Carmel  River.  Sand  is  introduced  by  the  river  primarily  during  the 
winter  and  spring  when  rainfall  is  sufficient  to  maintain  a  flow.  During 
the  summer  and  fall,  the  beach  closes  the  river  and  a  brackish  lagoon 
forms  behind  the  beach,  thus  ceasing  the  river-carried  sand  supply.  A 
secondary  source  of  sediment  is  the  shoreline -outcropping  rock. 

When  the  winter  rains  have  supplied  a  sufficient  amount  of  precipi- 
tation to  produce  a  steady  flow  in  the  river  and  the  lagoon  is  filling, 
the  mouth  of  the  river  is  opened  by  a  bulldozer  (Fig.  5>)«  Failure  to 
open  the  river  at  the  proper  time  could  result  in  flooding  of  homes  at 
the  lower  end  of  the  Carmel  Valley.  This  operation  is  performed  under 
the  direction  of  the  head  ranger  at  Point  Lobos  State  Reserve  and  normal- 
ly occurs  once  or  twice  each  winter.  During  the  winter  of  1 966-1 967* 
this  operation  was  repeated  2  5  times  because  of  alternating  rainy  and 
dry  periods.   The  possibility  exists  that  heavy  rains  coupled  with  heavy 
surf,  a  high  tide,  and  a  strong  onshore  wind  might  prevent  this  opening 
and  flooding  would  occur  [Mr.  Don  Rich,  personal  communication] . 

One  cut  is  all  that  is  required  as  the  flow  of  the  river  widens  the 
cut  until  the  river  is  at  least  100  ft  wide  at  the  beach  (Fig.  6).   The 
author  estimated  that,  at  the  point  of  maximum  opening,  10,000  to  20,000 
cubic  yard  of  sediment  had  been  removed  from  the  beach. 

Much  of  the  fine  sediment  material  was  observed  moving  out  into  the 
bay  and  was  obviously  lost  to  the  study  area.  The  coarse  sediment  un- 
doubtedly settled  quickly  and  was  either  moved  with  the  littoral  drift 
or  redeposited  on  the  beach  by  wave  action  when  the  river  mouth  was 
closed  again. 
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Figure  6.   Mouth  of  the  Carmel  River  After  Flow  Has  Widened 
the  Bulldozer  Cut. 


C.   PREVIOUS  WORK 

The  geology  of  the  area  in  the  vicinity  of  Carmel  Bay  has  been  de- 
scribed in  some  detail  by  Lawson  [ 1 893  ]  j  Beal  [l 91 5] ,   Bowen  [l 96^] j 
Nili-Esfahani  | 1 965 J  and  Simpson  [l 972] .   Of  interest  to  this  study  are 
the  formations  that  are  in  the  study  area  or  are  associated  with  the 
Carmel  River  and  San  Jose  Creek  waterbeds.  The  Santa  Lucia  porphyritic 
biotite  granodiorite  (Cretaceous)  is  apparently  the  most  important  source 
of  sediments  for  the  area.  The  Carmelo  Series  (Paleocene),  interbedded 
pebble  conglomerate,  sandstone,  and  siltstone,  and  Monterey  shale  (Upper 
Miocene),  a  silicious  shale,  provide  a  smaller  amount  of  sediments.  Sur 
Series  metamorphic  gneisses  (pre-Cretaceous )  apparently  contribute  sig- 
nificantly to  the  heavy  minerals  supplied  to  the  area  as  they  dominate 
the  upper  regions  of  the  Carmel  River  watershed  [California  State  Depart- 
ment of  Water  Resources,  I960]. 

Diver  description  of  the  study  area  has  been  almost  nonexistent. 
Numerous  dives  have  been  conducted  in  and  around  the  head  of  the  Carmel 
Submarine  Canyon.   Shepard  and  Dill  [1966],  McLean  and  Peckham  [1961], 
Moritz  [1968],  and  Wallin  [1968]  have  provided  descriptions  of  San  Jose 
Creek  beach  and  the  head  of  the  canyon.  Sand  falls  and  slump  scars  have 
been  noted  by  Moritz . 

Bathymetric  surveys  of  the  bay  were  accomplished  by  the  U.S.  Coast 
and  Geodetic  Survey  (USCGS)  in  1 883  (No.  l5U8a)  and  1933  (No.  K$k53)> 
Shepard  and  Emery  [19U1]  have  conducted  a  lead  line  survey  of  the  head 
of  the  Carmel  Submarine  Canyon.  Zardeskas  [1971]  conducted  the  first 
bathymetric  survey  using  electronic  navigation  and  echo-sounding  equip- 
ment.  This  survey  is  considered  to  be  more  accurate  than  previous  sur- 
veys but  does  not  provide  sufficient  detail  inside  the  1 0-f m  contour. 
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Studies  of  the  currents  in  the  bay  are  almost  non-existent.  Tides 
of  the  bay  are  the  mixed  type  that  is  characteristic  of  the  Pacific  Coast. 
Diurnal  differences  between  mean  higher  high  water  and  mean  lower  low 
water  average  5>«2  ft  [U.S.  Coast  and  Geodetic  Survey,  1  971 J  . 

Bascom  [1°6U]  conducted  a  study  of  the  seasonal  changes  in  the  pro- 
file of  the  Carmel  River  State  Beach.  These  seasonal  changes  were  of 
little  significance  in  this  study.  However,  these  seasonal  cycles  show 
that  sand  is  moved  from  the  beach  to  an  offshore  bar  by  winter  storms 
and  is  carried  back  onto  the  beach  by  summer  swell. 

Carter  [ 1971  J  studied  the  marine  sediments  of  the  bay  using  grain- 
size  analysis  techniques  and  Griffin  [1969J  has  conducted  a  heavy-mineral 
analysis  of  the  beaches  in  the  bay. 
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II.  METHODS  OF  DATA  ACCUMULATION 

A.   BATHYMETRY 
1 .   Survey 

A  bathymetric  survey  of  the  area  was  conducted  on  May  12,  1972. 
This  survey  was  undertaken  in  order  to  provide  information  for  compari- 
son with  previous  surveys  by  the  USCGS  and  the  Naval  Postgraduate  School 
(NPS).   Large  depth  discrepancies  between  the  various  surveys  could  pro- 
vide an  indication  of  slumping.  Additionally,  the  survey  provided  an 
accurate  nearshore  chart  for  wave  refraction  diagrams. 

Four  transit  locations  along  the  beach  in  the  study  area  labelled 
A,  B,  C,  and  D  (shown  on  the  finished  chart,  Fig.  13)  were  located  using 
the  U.S.  Geological  Survey  Horizontal  Control  Station  "Whaler's  Knoll 
Number  2"  (Fig.  2),  the  new  Carmel  Mission  Spire,  and  the  references 
used  in  the  1  933  USCGS  survey  that  were  designated  "Hudson"  and  "Granite 
Point"  on  the  finished  chart.   Two  Dietzgen  6000  Series  land  surveying 
transits  were  used,  first  at  locations  A  and  B  and  last  at  locations  C 
and  D,  to  determine  the  position  of  the  survey  boat.  A  large  marker  was 
placed  on  Whaler's  Knoll  Number  2  which  could  be  seen  from  the  transit 
locations.   This  served  as  a  zero  reference  for  the  transit  operations 
and  all  angles  were  measured  relative  to  this  marker. 

An  Apelco  Model  MR-201B  fathometer  was  mounted  on  the  transom  of 
a  13 -ft  small  craft  powered  by  a  9-5>-hp  outboard  motor.  An  aluminum 
pole  painted  international  orange  was  mounted  vertically  above  the  fath- 
ometer transducer.  Communication  between  the  survey  boat  and  the  transit 
operators  was  maintained  with  1 -w  two-way  radios  at  each  location.   The 
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position  of  the  survey  boat  was  determined  at  least  every  minute  during 
the  survey. 

Track  lines  were  intended  to  be  perpendicular  to  the  beach  at 
30-yard  intervals.   This  was  difficult  to  maintain  because  of  large 
amount  of  kelp  (Microcystis  pryifera)  and  difficulty  in  seeing  prelocated 
beach  reference  marks . 

2.  Depth  Corrections 

All  depths  on  the  finished  chart  were  corrected  to  mean  lower 
low  water  datum.  In  order  to  correct  for  the  depth  of  the  transducer, 
differences  in  calibrated  sound  velocity ,  and  instrument  error,  compari- 
sons between  the  fathometer  readings  and  lead-line  readings  were  made 
periodically  over  a  known  flat  bottom. 

A  tidal  correction  was  applied  using  the  NPS  tide  gage  located 
on  Monterey  Municipal  Wharf  Number  Two  in  Monterey  Bay  (just  north  of 
Carmel  Bay)  as  a  reference.  Tidal  differences  between  Monterey  Bay  and 
Carmel  Bay  rarely  exceed  0.£  ft  [Dr.  W.C.  Thompson,  personal  communica- 
tion]. Since  the  fathometer  sound  beam  spreads  spherically  within  its 
beam  width,  a  correction  is  required  if  the  bottom  is  not  horizontal. 

True  slope  of  the  bottom  was  obtained  by  plotting  the  various 
data  points  and  measuring  the  distance  between  two  points  across  the 
slope .   The  fathometer  readings  are  known  at  each  point  and  the  slope  is 
easily  obtained.  A  correction  for  slope  was  applied  using  tables  prepared 
by  Zardeskas  [ 1 971 ]  . 

B.   DIVER  OBSERVATIONS 

An  examination  of  the  study  area  was  made  over  a  7-month  period 
using  SCUBA,  during  which  23  dives  were  made  in  the  area.  A  log  of  six 
of  the  more  interesting  dives  is  contained  in  Appendix  A.  On  all  dives 
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an  effort  was  made  to  record  general  bottom  types,  sand  thicknesses, 
unusual  features,  bottom  slopes  and,  when  applicable,  sediment  and  rock 
samples  were  taken. 

C.   LITTORAL  DRIFT 

Littoral  drift  was  qualitatively  determined  using  statistical  wave 
data  from  National  Marine  Consultants  (NMC)  [i960],  and  wave  refraction 
methods  given  in  U.S.  Naval  Hydrographic  Office  [195>8]  Publication  Number 
H.O.  23U •  No  attempt  was  made  to  quantitatively  determine  the  longshore 
sand  transport  because  of  insufficient  wave  data  and  because  empirical 
relationships  to  determine  volume  transported  are  based  on  a  smooth  and 
straight  beach  [Bowen  and  Inman,  1 966 ] .   The  beach  in  the  study  area  is 
broken  in  numerous  places  by  large  granodiorite  outcrops. 

1  .   Statistical  Wave  Data 

The  wave  data  used  are  the  results  of  hindcasting  from  weather 
map  analyses  for  the  three  years  1  95>6,  1957  and  1 958  (fig-  7).  Sea  and 
swell  data  are  given  by  period,  deep  water  wave  height,  frequency  of 
occurrence  and  direction.  Direction  is  given  in  22.5°  increments 
(¥,  W1TW,  NW,  etc.). 

None  of  the  seven  stations  along  the  California  Coast  for  which 
the  data  are  presented  were  located  in  the  vicinity  of  Carmel  Bay.   The 
data  of  the  nearest  station  to  the  north  (Number  3),  located  approximate- 
ly 50  nm  west  of  the  entrance  to  San  Francisco  Bay  (Latitude  37.6  N, 
Longitude  123-5  W),  and  the  nearest  station  to  the  south  (Number  h) 
(Latitude  35-5°  W,  Longitude  122.0°  W),  showed  no  significant  differences. 
For  this  reason,  the  data  from  Station  h  were  used. 

The  NMC  information  is  incomplete  because  it  is  based  on  synoptic 
meteorological  data  from  the  North  Pacific  and  does  not  consider: 
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Figure  7-  Swell  Data  from  National  Marine  Consultants 
for  Station  k-    (after  Cherry,  1965). 
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a.  local  waves  generated  by  the  diurnal  sea  breezes,  and 

b.  southerly  swell. 

The  first,  local  waves,  could  be  of  considerable  significance  in 
sand  movement.   The  second,  because  of  the  great  amount  of  refraction 
that  would  be  involved,  would  probably  not  contribute  significantly  to 
the  sand  transport . 

The  data  show  that  over  b0%   of  the  deep  water  wave  energy  that 
arrived  off  the  coast  comes  from  the  directions  of  west,  west -northwest 
and  northwest.   Over  h.0%   of  the  total  energy  comes  from  the  northwest 
alone   [ Cherry,  1 96£ ] • 

2 .   Construction  of  Wave  Refraction  Diagrams 

Wave  refraction  diagrams  were  constructed  for  waves  from  the 
west,  west-northwest  and  northwest  by  the  author  and  by  previous  inves- 
tigators [Wallin,  1968  J.  The  methods  given  in  H.O.  23a  were  used.  All 
calculations  were  accomplished  using  a  computer  program  that  eliminates 
the  need  for  tables  to  determine  the  ratio  of  the  depth  of  the  water  to 
the  deep  water  wave  length.  A  copy  of  the  computer  program,  formats  for 
data  cards  and  a  sample  of  the  output  are  listed  in  Appendix  B. 

Direction  of  transport  for  San  Jose  Creek  Beach,  located  just 
south  of  the  study  area,  and  Carmel  River  State  Beach  were  determined. 
Refraction  diagrams  were  considered  more  accurate  than  previous  diagrams 
because  of  the  improved  bathymetric  data  provided  by  Zardeskas  [1971 J 
and  the  author's  own  nearshore  survey. 

A  wave  refraction  diagram  for  swell  from  the  west-northwest  with 
a  period  of  11  sec  is  shown  in  Fig.  8. 
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Figure  8.  Wave  Refraction  Diagram  for  West-Northwest  Swell  with 
a  period  of  1 1  sec . 
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D.   CURRENT  STUDIES 

1  .  Savonius  Rotor  Current  Meter 

Two  attempts  were  made  to  measure  bottom  currents  outside  the 
surf  zone.   The  first  involved  the  use  of  a  Hydro  Products  Model  501  in 
situ  current  recording  system.  This  instrument  measures  current  speed 
with  a  Savonius  rotor  and  current  direction  with  a  direction  vane.  Di- 
rection ,  speed,  and  temperature  are  recorded  with  a  Rustrak  recorder. 
Current  speed  in  shallow  water  can  be  measured  with  this  sytem  only  if 
the  water  surge  component  is  negligible  compared  to  the  current  [Dr.  E.B. 
Thornton,  personal  communication].   For  this  reason,  the  meter  was  operated 
for  1  week  in  the  spring  with  the  hope  that  1  or  2  days  of  minimal  wave 
activity  might  ensue . 

The  meter  required  a  vertical  orientation  in  order  to  eliminate 
any  gravity  component  on  the  directional  vane  when  recording.  Since  the 
bottom  slope  was  1  k°   at  the  chosen  location,  the  vane  was  fitted  with  a 
small  glass  floatation  chamber  to  eliminate  this  problem  when  submerged 
in  ■  salt  water.  Additionally,  the  meter  was  equipped  with  a  plexiglass 
cover  for  diver  observation  of  the  recorder. 

A  ££0-lb  rectangular  (1  8  inch  x  18  inch  x  2  inch)  lead  base  im- 
bedded with  stainless  steel  studs  for  meter  attachment  was  prepositioned 
using  the  R/V  ACANIA  in  60  ft  of  water  in  the  center  of  the  study  area. 
One  of  the  sediment  samples  (Number  6,  Fig.  11 )  was  taken  next  to  this 
base .   The  meter  was  lowered  from  a  small  boat  and  bolted  to  the  base  by 
divers  (Fig.  9). 

2 .   Ducted  Current  Meter 

The  second  effort  used  a  Marine  Advisor's  ducted  current  meter 
system  attached  to  a  base  fabricated  at  the  machine  facility  of  NPS 
(Fig.  10).  This  system,  unlike  the  Savonius  rotor  model,  can  essentially 
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Figure  9.  Divers  In      >g  Savon 


Figure  10.  Diver  Placing  Ducted  Current  Meter  and  Base  on  the  Bottom. 
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eliminate  the  wave  surge  component  simply  by  physically  aligning  the 
axis  of  the  meter  parallel  to  the  crests  of  the  wave-induced  sand  ripples. 
Its  major  disadvantage  is  the  requirement  for  surface  recorders.  The 
power  for  the  recorders  was  supplied  by  a  12-v  automobile  battery  and  a 
DC /AC  inverter.   The  recorders  consisted  of  a  ducted  current  meter  read- 
out Model  S-6A  (0-1 5>  knot  range)  and  a  Varian  Model  G-11A  strip  chart 
recorder  adjusted  for  a  full  scale  deflection  of  2  knots.   The  meter  was 
placed  in  60  ft  of  water  directly  seaward  of  the  Carmel  River  lagoon. 
The  axis  of  the  ducted  meter  was  oriented  so  that  the  wave  surge  component 
was  minimal.  Data  was  taken  for  periods  of  approximately  30  min. . 

E.   SEDIMENT  SIZE  ANALYSIS 

1  .   Collection  and  Analysis  Procedures 

Eighteen  sediment  samples  were  taken  in  the  area  and  grain-size 
analyses  were  conducted  using  the  procedures  outlined  in  Krumbein  and 
Pettijohn  [1938  J.  Twelve  of  the  samples  were  taken  by  the  author  using 
SCUBA,  five  were  taken  at  selected  locations  along  the  beach,  and  one 
(sample  Number  18)  was  taken  in  the  bed  of  the  Carmel  River  in  the  vi- 
cinity of  the  U.S.  Geological  Service  River  Gaging  Station  near  the 
bridge  of  California  Highway  1  (Fig.  11).  All  diver  samples,  with  the 
exception  of  sample  Number  12,  were  taken  across  the  crests  of  sand 
ripples  to  insure  uniformity  (Fig.  12).   Samples  Number  11  and  12  were 
taken  in  a  crest  and  adjacent  trough  respectively  for  comparison. 

From  each  sample  a  subsample  of  approximately  IjO  g  was  taken. 
The  subsample  was  washed  with  distilled  water  in  order  to  remove  the 
salt,  allowed  to  settle,  and  the  excess  water  was  decanted.  The  sub- 
sample  was  then  wet  sieved  through  a  stainless  steel  U0   screen.  No  fine 
fraction  (  >  U0)   of  any  significance  was  noted  for  any  of  the  samples. 
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Figure  11.   Locations  of  Sediment  Samples. 
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Figure  12.  Diver  Taking  Sediment  Sam] 


The  coarse  fraction  was  placed  in  a  beaker,  dried  and  weighed.  After 
weighing,  the  coarse  fraction  was  dry  sieved  through  8-inch  U.S.  stan- 
dard screens  at  0.^0  intervals  using  a  Ro-Tap  automatic  shaker.  The 
fraction  retained  on  each  of  the  screens  was  then  weighed  to  the  nearest 
0.1  mg.  Finally,  the  cumulative  weight  of  fractions  was  compared  to  the 
original  weight  to  check  for  gross  errors. 
2 .   Computer  Analysis 

Folk  and  Ward  [1  95>7j  statistical  parameters  were  calculated 
using  an  IBM  360  computer  and  a  computer  program  prepared  by  ¥.  R. 
Anikouchine  [Dinger,  1970,  p.  31  ]  and  slightly  modified  by  Miss  Sharon 
Raney  (NPS  Computer  Facility).  A  copy  of  the  computer  program,  formats 
for  data  cards  and  a  sample  of  the  output  are  listed  in  Appendix  B. 
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III.   RESULTS 

A.  BATHYMETRY 

Chart  preparation  was  accomplished  by  first  preparing  a  smooth  sheet. 
This  required  the  plotting  of  all  fixes  that  had  been  taken  on  a  1 :5000 
scale  sheet j  printing  the  corrected  depths  on  the  plotted  positions  and 
contouring  the  sheet  in  increments  of  10  ft.   The  contours  range  from 
20  ft  to  100  ft.   The  1 0-f t  contour  was  omitted  because  of  insufficient 
data.   Moderate  surf  activity  on  the  day  of  the  survey  precluded  boat 
operations  very  near  shore.  The  shoreline  for  the  chart  was  obtained 
from  the  USCGS  1933  smooth  sheet. 

The  finished  chart  (Fig.  13)  was  prepared  on  tracing  paper  from  the 
smooth  sheet.   It  has  been  reduced  for  inclusion  in  this  paper. 

Of  interest  is  the  shoreward  trend  of  the  contour  lines  between 
transit  locations  B  and  C.  This  might  possibly  be  due  to  a  deprivation 
of  sediment  supply  caused  by  a  large  granodiorite  outcrop.  Additionally, 
there  exists  a  slight  seaward  extension  of  contour  lines  just  offshore 
of  transit  location  C.   This  might  be  caused  by  a  larger  rate  of  sediment 
supply  and/or  an  underlying  formation.  This  seaward  bend  of  contour 
lines  is  dramatically  shown  in  Fig.  2. 

Comparison  of  the  finished  chart  with  the  USCGS  1883  and  1933  surveys 
revealed  no  significant  depth  changes.   The  study  areas  of  this  chart  and 
that  of  Shepard  and  Emery  [ 1 9Ul ]  did  not  sufficiently  overlap  and  thus 
provided  little  useful  correlation.  The  1971  survey  (Fig.  2)  did  not, 
because  of  limitations  of  the  towed  sounder,  survey  this  close  to  shore. 
The  10-fathom  curve  of  the  1971  survey  was  taken  from  the  USCGS  1933 
survey. 
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Figure  13-  Nearshore  Bathymetric  Chart, 
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Shepard  and  Emery  had  previously  compared  their  survey  and  the  USCGS 
1883  survey  and  shallower  soundings  near  the  canyon  head  (south  of  the 
study  area)  indicated  that  a  slide  or  a  slump  had  occurred. 

B.  DIVER  OBSERVATIONS 

The  area  of  study  was  broken  into  Subareas  1  through  7  (Fig-  1ii)  in 
order  to  describe  observations  by  the  author  using  SCUBA. 

Subarea  1  was  entirely  granodiorite  from  the  shoreline  to  a  depth  of 
65  ft.   There  appears  to  be  no  way  that  sediments  could  be  transported 
across  this  barrier.  Wave  refraction  studies  [Dr.  W.C.  Thompson,  per- 
sonal communication]  and  heavy  mineral  analysis  [Griffin,  1969]  reinforce 
this  assumption. 

Subarea  2  is  almost  entirely  covered  with  a  thick  layer  (greater  than 
5  ft)  of  fine  to  coarse  sand.   One  small  granodiorite  outcrop  exists  in 
approximately  25  ft  of  water  just  west  of  the  Carmel  River.  Bottom  slopes 
in  the  area  range  from  nearly  horizontal  in  shallow  areas  to  25°  -30°  in 
the  southwestern  corner  of  the  area. 

Subarea  3  is  a  combination  of  rock  and  sand  bottom.  There  is  no 
physical  barrier  to  sand  movement  in  any  direction  in  this  area. 

Subarea  i|.  is  entirely  sand  with  the  exception  of  one  large  granodio- 
rite outcrop.  This  outcrop  is  narrow  (25  ft  to  50  ft)  and  long  (approxi- 
mately 100  yards)  and  is  oriented  parallel  to  the  beach.   On  the  western 
side  of  this  unbroken  outcrop,  bottom  slopes  range  from  5°  to  20°.   The 
bottom  on  the  eastern  side  has  very  gentle  slopes.  Sand  sizes  on  the 
western  side  are  noticeably  finer  than  on  the  eastern  side.   This  large 
outcrop  precludes  any  significant  sediment  transport  in  a  westerly 
direction.  Sediment  movement  in  the  area  appears  to  be  confined  to  the 
area  between  the  beach  and  the  outcrop  and  is  generally  in  a  southerly  di- 
rection. The  grain  size  of  the  sand  gave  evidence  of  much  local  erosion. 
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Figure  U|.  Subareas  Surveyed  Using  SCUBA. 
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Subarea  J>  contains  the  southern  end  of  the  granodiorite  outcrop  de- 
scribed in  Subarea  i|.  The  outcrop  is  broken  in  several  locations  near 
the  end.   Through  these  breaks  there  was  evidence  of  sediment  movement 
as  large  fans  or  aprons  of  sedimentary  material  existed  at  the  base  of 
each  of  the  breaks.  Some  sediment  movement  was  observed  by  the  author. 
It  is  presumed  that  a  large  swell  is  required  before  large  amounts  of 
sediment  are  moved.  The  remainder  of  the  area  is  covered  with  a  very 
thick  layer  of  very  coarse  sand. 

Subarea  6  is  located  just  seaward  of  a  large  rock  outcrop  that  is 
west  of  the  entrance  on  San  Jose  Creek.  This  outcrop  is  awash  about  3 
to  k   ft  at  low  tide  and  is  called  the  "washrock"  by  local  divers.  This 
area  contains  the  head  of  the  Carmel  Canyon.  Slopes  in  this  area  are  no 
less  than  J4O0.   The  bottom,  although  sides  would  probably  be  a  more 
appropriate  term,  is  entirely  made  up  of  large  granodiorite  boulders. 
Depths  in  the  area  reach  1 20  ft  within  a  few  tens  of  yards  from  the 
"washrock" . 

Subarea  7  is  almost  entirely  rock  covered.   There  appears  to  be 
little  or  no  sediment  transport  across  this  area. 

Subarea  5>  is  the  most  interesting  area  observed.   Previous  studies 
have  described  "sand  chutes"  along  the  north  wall  of  the  canyon  head 
[Wallin,  1968].  The  author  noted  some  minor  changes  in  the  depths  of  the 
southern  end  of  the  subarea  over  the  several  dives  that  were  accomplish- 
ed. Significant  changes  in  depth  (greater  than  2$   ft)  have  been  noted 
by  Cdr.  Don  Ferrin,  USN  (Ret.)  [personal  communication  J  several  times 
over  the  past  12  years.   The  subarea  is  almost  entirely  covered  with 
coarse  sand.  Much  of  this  appears  to  have  been  derived  from  outcrops 
in  the  near  area . 
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C.  LITTORAL  DRIFT 

Wave  refraction  diagrams  constructed  indicated  that  the  littoral 
drift  in  the  study  area  is  south.  Additionally,  the  drift  along  San 
Jose  Creek  Beach  is  to  the  north  and  a  node  exists  just  north  of  the 
head  of  the  Carmel  Canyon.  The  area  of  this  node  was  described  in  a 
previous  section  (Subarea  £) • 

Stream  mouths  generally  tend  to  migrate  in  the  direction  of  littoral 
drift  [U.S.  Army  Coastal  Engineering  Research  Center,  1 966] .  Evidence  of 
the  southerly  drift  in  the  study  area  was  given  by  a  southerly  migration 
of  the  mouth  of  the  Carmel  River  after  opening.   The  rate  of  this  drift 
was  observed  to  be  dependent  upon  the  wave  activity  and  angle  of  inci- 
dence of  the  waves  upon  the  beach.  Further  southward  migration  of  the 
river  mouth  is  prevented  by  rock  outcrops.   Evidence  of  the  northerly 
drift  along  San  Jose  Creek  Beach  is  given  by  the  long-term  northward 
migration  of  the  mouth  of  San  Jose  Creek  (Fig.  3). 

D.  CURRENT  STUDIES 

Current  data  obtained  from  the  Savonius  rotor  meter  was  judged  to  be 
invalid.  Moderate  to  heavy  swell  persisted  during  the  entire  week  that 
the  meter  was  in  place.   The  bottom  wave  surge  induced  by  this  swell  pro- 
duced an  excessive  amount  of  scatter  and  an  unreadable  record  resulted. 

The  ducted  meter  produced  better  results  than  the  Savonius  rotor 
model.   Currents  at  the  chosen  station  for  the  period  that  the  meter  was 
on  the  bottom  were  less  than  0.2  knots  and  this,  because  of  its  oscilla- 
tory nature,  was  judged  to  be  a  component  of  wave  surge.   Currents  in 
water  less  than  20  ft  depth  could  not  be  obtained  because  of  the  hazards 
to  divers  and  equipment  in  the  surf  zone. 
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Even  minimal  currents  can  be  detected  by  a  diver  who  is  required  to 
swim  against  them.  During  the  entire  study  with  the  more  than  20  dives 
that  were  completed,  no  currents  were  noted  except  when  entering  and 
leaving  the  surf  zone.   This  observation  reinforces  the  results  obtained 
using  the  ducted  meter. 

E.   SEDIMENT  SIZE  ANALYSIS 

The  tabulated  results  as  obtained  from  the  computer  program  are 
presented  in  Table  I. 
1  .  Mean  Grain  Size 

The  mean  grain  size  of  the  samples  taken  ranged  from  2.710  to 
-1  .580.  Six  of  the  samples  were  classified  as  fine  sand,  two  as  medium 
sand,  three  as  coarse  sand  and  seven  as  very  coarse  sand. 

Sand-silt-clay-gravel  relationships  were  computed  and  plotted 
on  a  sand-silt-gravel  diagram  similar  to  the  tetrahedron  scheme  of 
Krumbein  and  Sloss  [1963]  (Fig.  15)-   "Gravel"  refers  to  particles  with 
mean  grain  sizes  coarser  than  -1 .00  and  finer  than  -6.00. 

In  general,  mean  grain  size  was  larger  nearshore  and  in  the 
southern  end  of  the  study  area  (Fig.  16).  Data  from  Carter  [ 1971 J  were 
used  to  assist  in  the  locating  of  the  30  contour.  The  seaward  bend  of 
the  mean  grain  size  contours  in  Subarea  5>  probably  indicates  than  sedi- 
ments are  being  transported  seaward  in  this  area. 
2 .   Standard  Deviation 

The  values  of  standard  deviation  ranged  from  0 .I4I  0  to  1  .350. 
One  sample  was  classed  as  poorly  sorted.   Two  samples  were  well  sorted, 
and  1^  samples  were  moderately  sorted.  No  trends  were  noted  in  this 
statistical  measure.  A  mineral  analysis  of  the  samples  might  provide 
an  insight  into  this  parameter. 
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Figure  16.   Sediment  Mean  Grain  Size  Distribution, 
partially  on  data  from  Carter,  1 971 .  ) 
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3 .  Skewness 

Skewness  values  ranged  from  -O.I4.3  to  O.63.  Three  samples  were 
classified  as  very  negatively  skewed,  six  as  negatively  skewed,  five  as 
nearly  symmetrical  (one  of  these  was  sample  Number  1 8  that  was  taken  in 
the  bed  of  the  Carmel  River),  two  were  positively  skewed  and  two  very 
positively  skewed.  Sample  Number  11,  taken  across  a  ripple  crest,  was 
very  negatively  skewed  (-0.37).  Sample  Number  12,  taken  in  the  trough 
of  the  ripple  from  which  sample  Number  11  was  taken,  was  positively  skewed, 
This  can  be  readily  understood  since  positive  skewness  implies  skewed 
toward  the  fine  (more  coarse  material  than  fine)  and  negative  skewness 
implies  toward  the  coarse.   The  wave  surge  turbulence  is  greater  at  the 
crests,  thus  producing  negative  skewness. 
h>     Kurtosis 

Kurtosis  values  ranged  from  0.91  to  3«21.  One  sample  was  classi- 
fied as  extremely  leptokurtic,  four  as  very  leptokurtic,  three  as  lepto- 
kurtic  and  eleven  as  mesokurtic .   The  samples  taken  in  the  southern  half 
of  the  study  area  and  either  in  shallow  water  or  on  the  beach  face  tended 
to  be  of  a  leptokurtic  nature. 
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IV.   SUMMARY  AND  CONCLUSIONS 

Diver  observations  in  the  area  showed  that  the  sand  transport  in  the 
study  area  is  generally  southward.  This  was  both  observed  in  the  surf 
zone  and  surmised  from  the  evidence  of  previous  sand  movement  in  the 
southern  end  of  Subarea  h   and  the  northern  end  of  Subarea  £. 

The  bathymetric  survey  provided  an  indication  of  sand  deprivation  in 
Subarea  l\   and  a  slight  seaward  bending  of  depth  contours  in  Subarea  £. 
Comparisons  of  the  author's  bathymetric  chart  with  previous  charts  pro- 
vided no  real  evidence  of  slumping  or  sliding.  Previous  studies  have 
given  some  evidence  for  this  [Shepard  and  Emery,  19hl] • 

Littoral  drift  was  determined  by  wave  refraction  diagrams  to  be  south 
in  the  study  area  and  north  along  San  Jose  Creek.  A  nodal  point  for  sand 
transport  exists  just  north  of  the  head  of  the  Carmel  Canyon  (Fig.  1?)- 

Sediment  size  analysis  reinforced  the  existence  of  a  nodal  point  for 
sediment  transport. 

The  two  major  sources  of  sedimentary  material  for  the  study  area  ap- 
pear to  be  the  Carmel  River  and  local  erosion  of  granodiorite  outcrops. 
The  Carmel  River  provides  sediments  only  during  the  winter  when  precipi- 
tation is  sufficient  to  create  a  flow.   Local  erosion  of  outcrops  is 
continuous  but  larger  during  the  winter  when  winter  storms  produce  in- 
creased wave  activity. 

Materials  provided  by  the  Carmel  River  are  moved  south  with  the 
littoral  drift,  mixed  with  locally  derived  sediments  and  deposited  along 
the  north  wall  of  the  Carmel  Canyon.  Sediments  provided  by  local  erosion 
of  outcrops  south  of  San  Jose  Creek  Beach  are  moved  north  with  the 
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Figure  1?.  Summary  of  Sand  Transport  in  the  Area. 
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littoral  drift,  irrbced  with  the  small  amount  of  sediments  provided  by 
San  Jose  Creek,  and  injected  into  the  Carmel  Canyon. 

Winter  storms  probably  provide  the  necessary  energy  to  induce  gravity 
sliding  or  slumping  of  deposited  material.   This  material  is  dumped  into 
the  head  of  the  Carmel  Submarine  Canyon  for  eventual  deposition  in  the 
deep  ocean. 
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V.   SUGGESTIONS  FOR  FURTHER  STUDIES 

The  following  studies  are  presently  being  conducted  within  the  bay: 
1  .   sediment  migration  within  Whaler's  Cove  (L.  Leopold,  San  Jose 

State  College,  in  progress); 
2.   gravity  survey  of  Carmel  Bay  (A.  Souto,  NPS,  in  progress). 
Future  beneficial  studies  of  the  marine  environment  of  the  bay  might 
include : 

1 .  currents  and  tagged  sand  studies; 

2.  seismic,  refraction  measurements; 

3.  heavy  mineral  analysis; 

1±.  carbon,  carbonate,  and  organic  nitrogen  analysis  of  sediments; 
5.  a  sand  budget  for  the  bay. 
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APPENDIX  A:   SELECTED  DIVER  LOGS 

DIVERS:  Howell,  Cronyn  DATE:   16  October  1971 

Subareas  h   and  5>  were  examined  during  this  dive. 

The  divers  entered  the  bay  through  the  surf  via  a  pocket  beach  ap- 
proximately 300  yards  south  of  the  Carmel  sewer  outfall.  Approximately 
50  yards  from  the  beach  is  an  outcrop  of  granodiorite  that  is  about  20 
yards  wide  and  ranges  from  1 0  to  1|0  ft  in  depth.   It  is  heavily  encrusted 
with  numerous  encrusting  organisms.   The  outcrop  is  lengthy  and  unbroken 
in  its  north-south  direction  and  appears  to  act  as  a  barrier  to  the  sand 
not  allowing  beach  sand  to  be  transported  seaward.  Seaward  of  the  outcrop 
the  sand  is  noticeably  finer  and  slopes  gently  downward  until  a  depth  of 
60  ft,  at  which  depth  there  was  a  marked  increase  in  the  slope . 

DIVERS:  Howell,  Mllward  DATE:  23  October  1971 

Subarea  3  and  the  southern  end  of  Subarea  2  were  examined  during  this 
dive. 

The  divers  entered  the  bay  through  the  surf  via  a  pocket  beach  ap- 
proximately 300  yards  south  of  the  Carmel  sewer  outfall.   The  break  in 
the  slope  was  located  and  the  divers  followed  this  slope  in  a  northerly 
direction  for  about  200  yards.  A  5-ft  pole  was  pushed  into  the  sand  at 
regular  intervals  in  order  to  determine  sand  thickness.   In  all  cases 
the  sand  cover  was  greater  than  £  ft.  The  divers  then  surfaced  and  swam 
to  end  of  the  temporary  trestle  that  was  being  used  to  lay  the  Carmel 
sewer  outfall  line.   Seaward  of  the  trestle  the  bottom  was  similar  to 
other  sand  bottoms  in  the  area  with  small  ripple  marks  and  gentle  slopes 
with  an  increase  in  the  slope  at  approximately  60  ft  depth. 


DIVERS:  Howell,  Cronyn  DATE:  29  October  1971 

Subareas  6  and  7  were  examined  during  this  dive. 

The  divers  entered  the  bay  through  the  surf  in  the  vicinity  of  San 
Jose  Creek.  After  submerging  just  south  of  the  wash  rock,  a  very  steep 
slope  in  the  sand  was  noted.  Although  no  angle  measuring  device  was 
used,  the  slope  appeared  to  be  greater  than  25>  •  Minor  slides  were 
readily  started  with  just  a  kick  of  a  flipper.   The  western  portion  of 
the  wash  rock  is  heavy  granodiorite  boulders  with  a  very  steep  slope . 
No  leveling  out  was  noted  up  to  a  depth  of  120  ft. 

DIVERS:  Howell,  Cronyn,  Millward  DATE:  30  October  1971 

Subarea  6  was  examined  during  this  dive. 

The  divers  entered  the  bay  through  the  surf  in  the  vicinity  of  San 
Jose  Creek.   The  bottom  from  the  beach  to  the  wash  rock  was  checked  for 
the  sand  channels  that  were  previously  noted  by  Wallin  [1968J.   These 
channels  were  not  found  nor  were  any  significant  sand  deposits  noted 
directly  east  of  the  wash  rock.  The  dive  commenced  at  the  same  location 
as  the  dive  of  29  October  1971 •  A  thorough  investigation  of  the  western 
side  of  the  outcrop  revealed  no  possible  passages  for  sand  to  be  trans- 
ported.  The  sand  on  both  sides  of  the  outcrop  contained  large  amounts 
of  detritus  (crab  shell,  broken  kelp,  etc.)  and  supported  a  substantial 
tubeworm  population  (approximate  density  20  per  ft  ).  It  was  concluded 
that  most  of  the  sand  introduced  into  the  bay  must  be  lost  into  the 
canyon  north  of  San  Jose  Creek. 

DIVERS:  Howell,  Mellor  DATE:   30  November  1971 

The  divers  attempted  to  launch  the  boat  in  Whaler's  Cove,  Point  Lobos 
State  Reserve,  but  were  prevented  because  of  high  swell.   During  the 
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night  and  early  morning  these  swells  had  destroyed  approximately  one 
half  of  the  trestle  that  had  been  erected  for  construction  of  the  600-ft 
Carmel  sewer  outfall.   Large  planks (approximately  6  inches  x  8  inches) 
that  had  been  bolted  to  piles  driven  into  granite  were  destroyed.  Need- 
less to  say,  a  dive  was  considered  inappropriate.  It  was  noted  that  the 
high  tide  for  that  day  was  an  exceptionally  high  spring  tide.   This  in 
conjunction  with  the  large  swell  was  the  reason  for  the  destruction. 

DIVERS:  Howell,  Mellor  DATE:   9  December  1971 

Subarea  k  was  examined  during  this  dive . 

The  divers  entered  the  bay  via  the  boat  ramp  at  Whaler's  Cove  at 
Point  Lobos  State  Reserve  in  the  boat.  The  purpose  of  the  dive  was  to 
locate  a  position  for  placement  of  current  meters.  A  position  was  lo- 
cated in  65  ft  of  water  on  a  sand  bottom  with  a  slope  of  17°3  approxi- 
mately 300  yards  south  of  the  Carmel  sewer  outfall  pier.   The  bottom  was 
very  clean  with  fine  white  sand  and  small  wave  induced  ripples. 
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APPENDIX  B:   COMPUTER  PROGRAMS 

COMPUTER  DATA  CARD  FORMATS 

i 
WAVE  REFRACTION  PROGRAM 

Card  1  :  Scale  Card. 

Col.  1-10     chart  scale   (FORTRAN  F1 0.Jj.  format) 

Example : 

^833- 33  ~~ \ 

scale:  833-33  ft  =  1  inch 
Card  2:   Period  and  Crest  Interval 
Col.  1-2      wave  period  in  sec 
8-1 1  .    crest  interval 
Example : 

(vo         2~7o  \ 

period:   10  sec 
crest  interval:   2.0 

SEDIMENT  SIZE  ANALYSIS  PROGRAM 

Cards  1  &  2:   Title  cards 

Col.  1-80  of  each  card  contain  alphanumeric  information  to 
appear  at  the  top  of  the  output. 

Card  3 :       Identifier  for  sample 

Col.  1-9  cruise  number 
1 0-1 2  sample  number 

13-18  sample  type        Columns  1-28  can  contain  any 
19-20  depth  legal  keypunch  character. 

21 -22  month 
23-21;  day 
25-28  year 
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Col.  29 -3k   latitude   (XX  XX. XX) 

35-1*0  longitude   (XX  XX.XX)  (First  digit  of  longitude  is 

determined  from  octant.) 
^7-61  depth  from  top  of  core(XXXXX-) 
62-66  length  of  core  (XXXXX. ) 
79  octant  (see  below) 


Example : 


f 72-HOW-01    10DIVER  50  525197236317821  ^9 


^ 


Cruise  Number: 

72-HOW-01 

Sample  Number: 

10 

Sample  Type: 

DITOR 

■ .  Depth : 

50 

Date: 

5/25/1  972 

Latitude : 

36°31 .78' 

Longitude : 

1 21 °55- 59' 

Card  h: 

Sample  Detail  Cards 

Col.  1*1  -I4.lt  phi  size  (absolute  value)  Fl*.2  (decimal  assumed) 
k5     sign  of  phi  size  (+  or  -) 
50-56  fraction  weight  F7 • U  (decimal  assumed) 
80  end  of  data  flag 

=8  if  last  phi  size  for  this  sample 
=9  if  last  phi  size  for  all  samples 


r 


Example : 


1 00+  39U56 

Phi  Size:  +1  .0 

Sample  Weight:  3-9ll56  g 

Flag:  8 


"*) 
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Coding  for  Octant  of  Geographic  Position 


Longitude  East 


0 


100°E 


180°W 


Longitude  West 
100°¥ 


6 

$ 

1 

2 

8 

7 

3 

h 

90°N 


Latitude 
North 


0U 

Latitude 
South 


90°S 
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Sample  of  Computer  Output  for  Wave  Refraction  Program. 


WAVE    PERIOD    =     11    SECONDS 
N    =    1.00 


WAVE    LENGTH    =       619.52    FEET 
SCALE:    1     IN.    =    833.3    FEET 


DEPTH 

FM. 

FT. 

1 

6 

2 

12 

3 

18 

4 

24 

5 

30 

6 

36 

7 

42 

8 

48 

9 

54 

10 

60 

15 

90 

20 

120 

25 

150 

30 

180 

40 

240 

50 

300 

60 

360 

D/LO 


L/LO 


WAVE    LENGTH 


NL 


0.010 

0.244 

0.019 

0.341 

0.029 

0.413 

0.039 

0.472 

0.048 

0.522 

0.058 

0.566 

0.068 

0.604 

6.077 

0.639 

0.087 

0.670 

0.097 

0.698 

0.145 

0.807 

0.194 

0.878 

0.242 

0.924 

0.291 

0.953 

0.387 

0.981 

0.484 

0.990 

0.581 

0.993 

ACTUAL 
FT. 

CHART 
IN. 

CHART 

IN. 

150.86 

0.18 

0.18 

211.15 

0.25 

0.25 

255.90 

0.31 

0.31 

292.37 

0.35 

0.35 

323.38 

0.39 

0.39 

350.42 

0.42 

0.42 

374.36 

0.45 

0.45 

395.79 

0.47 

0.47 

415.12 

0.50 

0.50 

432.65 

0.52 

0.52 

500.00 

0.60 

0.60 

543.78 

0.65 

0.65 

572.17 

0.69 

0.69 

590.13 

0.71 

0.71 

607.69 

0.73 

0.73 

613.61 

0.74 

0.74 

615.46 

0.74 

0.74 
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Sample  of  Computer  Output  of  the  Sediment  Size  Analysis  Program. 


MARINE    SEDIMENTS    NEAR    CARMEL    RIVER    STATE    BEACH,    CARMEL,     CALIFORNIA 

BUFORD    F.     HOWELL 

CRUISE    72-H0W-01    SAMPLE    NUMBER  5    DEPTH    20 

SAMPLER    TYPE    DIVER       DATE       5/25/1972       LAT.     36-32. 09N         LONG.    121-55. 72W 
DEPTH    FROM   TOP    OF    CORE  0.    MM.  LENGTH    CF    CORE  0.     MM. 

PHI  SAMPLE  FRACTION  ACCUMULATED 

SIZE  WEIGHT  PERCENT  PERCENT 

-3.00  1.9983             4.33  4.33 

-2.50  0.4966            1.08  5.41 

-2.00  1.3831             3.00  8.41 

-1.50  3.0940             6.71  15.12 

-1.00  10.5151  22.81  37.93 

-0.50  11.0373  23.94  61.87 

0.0  4.7591  10.32  72.19 

0.50  3.2197             6.98  79.18 

1.00  2.0965            4.55  83.72 

1.50  2.5737            5.58  89.31 

2.00  2.0364             4.42  93.72 

2.50  1.7354            3.87  97.60 

3.00  0.8432              1.83  99.43 

3.50  0.2150            0.47  99.89 

4.00  0.0499             0.11  100.00 

POST-ANALYTICAL     WEIGHT    IS       46.1032 

PHI    SIZES  AT    PERCENT    LEVEL  OF 

1     (C)  5  16  25  50     (M)  75  34  95 

-5.21  -2.68  -1.46  -1.21  -0.81  0.17  1.03  2.15    4    PT . 

-5.21  -2.68  -1.47  -1.25  -0.75  0.19  1.02  2.13    LINEAR 

SAND,     SILT,    CLAY,     RELATIONSHIPS 

GRAVEL  SAND  SILT  CLAY  TOTAL  SAND/MUD  TYPE 

37.93  62.07  0.0  0.0  100.00      9999.99  SAND 

TRASK    VALUES  4    PT .     PHIS 

Ql  Q2  Q3  SO  LOG    SO  SKG 

2.305  1.756  0.887  1.613  0.208  0.81 

INMAN    VALUES 

MEDIAN  MEAN  DEV.  SKEW.  2ND    SKEW.  KURT. 

-0.81  -0.22  1.24  0.48  0.44  0.94 

FOLK    AND    WARD    VALUE 

MEAN  -0.42  VERY    COARSE    SAND 

DEVIATION  1.35  POORLY     SORTED 

SKEWNESS  0.35  VERY    POSITIVELY    SKEWED 

KURTOSIS  1.44  LEPTGKURTIC 

TRASK    VALUES 
SI  Q2 

2.384  1.679 

INMAN    VALUES 

MEDIAN  MEAN 

-0.75  -0.23 

FOLK    AND    WARD    VALUE 
MEAN  -0.40 

DEVIATION  1.35 

SKEWNESS  0.31 

KURTOSIS  1.37 


LINEAR    1 
03 
0.876 

'HIS 

SO 
1.  649 

LOG    SO               SKG 
0.217             0.86 

DEV. 
1.25 

SKEW. 
0.42 

2ND    SKEW. 
0.38 

KURT 
0.93 

E 

VERY    COARSE    SAND 

POORLY    SORTED 

VERY    POSITIVELY    SKEWED 

LEPTOKURTIC 
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COMPUTER  PROGRAMS 

c 

C  THIS  PROGRAM  IS  DESIGNED  TO  OUTPUT  THE  NECESSARY 

C  INFORMATION  REQUIRED  FOR  MANUAL  CONSTRUCTION  OF 

C  WAVE  REFRACTICN  DIAGRAMS  USING  THE  METHOD  GIVEN  IN 

C  H.  0.  PUBLICATION  234.  THE  ONLY  INPUTS  REQUIRED 

C  ARE: 

C 

C  1  .  CHART  SCALE 

C  2.  WAVE  PERIOD 

C  3.  CREST  INTERVAL 

C 

C  THE    RATIO    L/LO    IS    PROVIDED    BY    A    SUBROUTINE    THAT 

C  USES    THE    NEWTON-S    ITERATION    METHOD    TO    SOLVE    THE 

C  EQUATION    L/L0=TANF(2*PI*D/L). 

C 

C 

CCMMGN/ IN  PIT/PERIOD, DEPTH 
READ(5,50)  SCALE 
■  50  FCPMATCF10.4) 
DO  1000  1=1,5 
READ(5,100) IPER,ZN 
100  FORMAT( I2,5X,F4.2) 
PERIOD=FLOAT( I  PER) 
C 

COMMENT:  CALCULATE  THE  DEEP  WATER  WAVE  LENGTH  (ZLO). 
C 

ZLC=5.12*IPER**2 
C 

COMMENT:  PRINT  A  TABLE  HEADING. 
C 

WR  ITE(6,200 ) I, IPERt ZLO tZN, SCALE 
200  FORMATt • 1 ' ///42X « TABLE ' ,I2,//15X' WAVE  PERIOD  =  «I2, 

*■  SECONDS' ,  lOXt  'WAVE  LENGTH  =  'F7.2,'  FEET'//15X'N  =  • 
*F4.2,26X,' SCALE:  1  IN.  =!,F6.1t'  FEET'/) 
WRITE(6,205) 
205  FCRMATC  •/,  18X,  '  DEPTH  «.8X,  •  D/LO *,  6X,  ■L/L  0 «  ,9X.  ■  WAVE  •  , 
*■  LENGTH"  ,7X, 'NL«  ,//  ,15X,C  FM.  «, 3X ,• FT .«, 27X, ' ACTUAL « 
*,'    CHART     CHART»/54X,«  FT.      IN.       IN.') 
K=0 
C 

COMMENT:  PROVIDE  DEPTHS. 
C 

DO  800  J=l,200 
210  IFtK.LE.9)  GO  TO  3CC 
220  IFIK.LE.25)  GO  TO  320 
230  IF(K.LE.SO)  GO  TO  330 
240  IF(K.LE.IOOO)  GO  TO  340 
300  K=K+1 

GO  TO  4CC 
320  K  =  K+5 

GO  TO  400 
330  K=K+10 

GO  TO  400 
340  K=K+25 
400  IFT=K*6 

RAT=IFT/ZLO 
C 

CCMMENT:IF  D/LO  IS  GREATER  THAN  .6,  GO  TO  THE  NEXT  TABLE. 
C 

IF    (RAT.GT..6)     GC    TO    1000 
CEPTH=FLOAT(  IFT  ) 
C 

COMMENT:CALL  SUBRCUTINE  TO  CALCULATE  L. 
C 

CALL    WAVE     (WAVED 
C 

COMMENT:    CALCULATE    L/LO. 
C 

RATIO=WAVEL/ZLO 
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c 

CCMMENT:  CALCULATE  DATA  FOR  CONSTRUCTION  OF  DIAGRAM. 
C 

CHART=WAVEL/SCALE 
ANL=ZN*CHART 
C 

CCMMENT:  PRINT  RESULTS. 
C 

60  0  hR I TE ( 6, 70 0)K, I  FT,  RAT, RATIO, WAV  EL, CHART, ANL 

700  FORMAT ( «0«,14X,I3,3X,I4,6X,F5.3,5X,F5.3,6X,F7.2,3X, 

*F4.2,5X,F4.2) 
800  CONTINUE 

WRITE(6,900) 
9C0  FCPMATC  '  ,  /  ,  130  ( 1H-)  ) 
10C0  CONTINUE 
STOP 
END 

SUBROUTINE  WAVE(WAVEL) 

CCN'MCN/INPUT/PERIOCCEPTH 

EXTERNAL  FCT 

DDPTH=DEPTH 

CALL  ZERO( WAVEL,FV,DERF,FCT,DDPTH,.005,20C,IER) 

RETURN 

END 

SUBROUTINE  Z ERO( X, F , DERF  ,FCT, XST , EPS , I  END  ,  I ER ) 
C 

C         PURPOSE 

COMMENT:     SOLVE    THE    GENERAL    NGN-LINEAR    EQUATIONS    CF    THE    FORM 
C  F(X)=0    BY    MEANS    OF    THE    NEWTCN-S     ITERATION    METHOD. 


C 


IER=0 
X  =  XST 
TOL  =  X 

CALL  FCT(TOL,F,DERF) 
TCLF=100.*EPS 
DO  6  I=1,IEND 
IF1F.E0.0.)  GO  TO  7 
IFtDERF.EQ.O.)  GO  TO  8 
2  CX=F/DERF 
X=X-DX 
TCL  =  X 

CALL  FCT(TOL,F,DERF) 
TCL=EPS 
A=ABS{X) 

IFIA.GT.1.)  TOL=TOL*A 
IF(ABS( DX) .GT.TOL )  GO  TO  6 
IF (ABS(F).LE.TOLF)  GO  TO  7 

6  CONTINUE 
IER=1 

7  RETURN 

8  IER=2 
RETURN 
END 


SUBROUTINE  FCT ( X , F , CER F ) 
COMMON /INPUT/PERI  CD,  DEPTH 

CATA  GG/32./,TW0PI/6. 28318/ 

ARG=DEPTH*TWOP I/X 

F=  GG*PERI0D**2/TWCPI-X/TANH(ARG) 

CERF  =  (-TANH(ARG)-ARG/COSH(  AR  G  )  *^:2  ) /TANH  (  ARG)**2 

RETURN 

END 


$k 


c 

C  SEDIMENT    SIZE    ANALYSIS    MAIN    PROGRAM 

C 

c 

COMMON    /BLK1/PHIS{8,2)  , N INM , NFAW , N SSC , SUMNL  ,NTRSK 

REAL    NEC/'-  «/ 

DATA  ANORTH/«N    «/,SOUTH/'S    '/iEAST/'E    '/, 
*WEST/'W    '/ 

DIMENSION  FRWT(IOO)  ,  PRCT ( 1 00  )  t CRMC (20 ) , STMC ( 20  ) , 
*EXMC(2  0) ,CPMB(20),STMB(2  0)  ,EXMB(20) , I CR ( 20  )  ,  JCR ( 20)  , 
*TATLE(20),TITLE{20) 

REAL* 8     CRUZ , CRUZ R,CRMB,CRMCt SMPLR 

DIMENSION  T(100)  ,PHI(  ICO) ,ACPC( 100) 

DATA    KKK/O/,MAT/O/,MBT/0/,MED/O/,KSM/0/ 

NIN0T=1 

902  N2=l 
NA=1 
NZ  =  1 
NL  =  1 
KB  =  1 
M=0 
PAWT=0. 

DO  903  K  =  l,.100 
T(K)=4.09 
ACPC(K)=999.99 
FRWT(K) =0. 0 

903  PHI(K)=99.99 
MA=0 

M8  =  0 

MBT=0 

PC=0 

MED  =  0 

SUMWT=0.0 

SUMPC=0.0 

ACPC( 1) =0.0 

READ(5,20)  TITLE 

20  FORMAT (20A4) 
READ(5,21)  TATLE 

21  FORMAT (20A4) 

1  NK  =  1 

2  READ(5,26)CPUZR,MCR,STATR,SMPLR,EXID,M0,DA,YR,LATA? 
*CEGLT, LNGA,DEGLNtPHIR,SIGN,BFWT,DFPTHR,CPLN»IQUD,NE 

26    F0RMAT(A7,A2,A3,A6T3A2  tA4  , 1  2  ,  F4.2  ,  12  ,  2F4.  2  ,  Al ,  4X, F7.4, 
*2F5.0, 12X, 211) 
TIMES  =  1  . 

IF(SIGM.EQ-NEG)     TIMES=-1. 
PHIR=(TIMES*PHIR) 

3  IF(PHIR.LT.3.9)  GO  TO  29 
IF(NZ.EQ.l)  GO  TC  48 
IF(PHIR.LE.4.1)  GO  TO  46 
NZ=2 

GO  TO  29 
46  NZ=1 

GO  TO  44 
48  IF(NL.FQ.2)  GO  TO  44 

NL=2 

NZ=2 
29  FRWTR=BFWT 

IF(FRWTR.GE„0. )  GO  TO  52 

IF( (FRWTR+. CD .GE.O.)  GC  TC  25 

WRITF{6,51)     FRWTR 
25    FRWTR=0.0 
52    PAI«\T  =  PAWT+FF.WTR 

SUMWT=PAWT 

GC    TO    49 
44    IF(NK.EQ.2)     GO     TO    59 

NK  =  2 

FWT1=BFWT*50. 

GO    TO    2 
59    FWT2=BFWT*50. 


60 
51 


53 


54 


FRWTR=FWT1-FWT2 

IFtFRWTR.GE.O. )    GO    TO    53 

IF( (FRWTR+O.Ol  }  .GE.O. J     GO    TO    55 

WRITE(6,51)     FRWTR 

FORMAT(25X,   'WEIGHING  ERROR.  FRACTION 
*F10.4,«  BEFORE  BEING  SET  TO  ZERO') 
c5  FRWTR=0.0 

FWT1=FWT2 

N I  NOT  =  2 

IF(NE.LT.8)  GO  TG  52 

PAWT=PAWT+FRWTR+FWT2 

SUMWT=PAWT 

K  =  K+1 

M  =  M+1 

KK=KK+1 

PHI(K)=PHIR 

FRWT(K)=FRWTR 

K  =  K+1 

M=M+1 

KK=KK+1 

FRWT(K)=FWT2 

PHI(K)=12. 

IFtPHIR.LE.il. 

GO  TO  100 
49  IF(NA.EQ.2).  GO 

NA=2 


WEIGHT  WAS  • 


)  PHI(K)=PHIR+1 
TO  50 


DETERMINE  OCTANT  FROM  FIRST  CARD  CF  SAMPLE 


DK=WEST 
IDH=LNGA 

IF(IQUD-2)  8,9,11 
IDH=IDH+100 
DG=ANORTH 
GO  TO  23 

IFUQUD-4)  12,13,15 
ICH=IDH+100 
DG=SOUTH 
GO  TO  23 
DK=EAST 

IF(IQUD-6)  17,16,19 
IDH=IDH+100 
CC=ANORTH 
GO  TO  23 
IF{ IQUD.GE 
DG=SOUTH 
CRUZ=CRUZR 
NCR=MCR 
STAT=STATR 
EXC=EXID 
DPTH=DEPTHR 
WRITE  (6,799)  TITLE 
FCRMATCl'  ,26X,20A4) 
WRITC(6,  800)  TATLE 
FORMAT ( »0« ,26X ,20A4//) 

WRITE  (6,831)  CRUZ, NCR, STAT, EXCSMPLR, MO, DA, YR, LATA, 
*DEGLT,DG, I DH ,DEGLN , DK ,DE PTHR ,CRLM 
801  F0RMAT(/36X, 'CRU ISE  ',A7,A2,'  SAMPLE  NUMBER 

SAMPLER  TYPE 


8 
9 

11 
12 
13 

15 

17 
16 

19 

23 


799 
800 


8)  IDH=IDH+100 


*•  DEPTH  «,A2//25X, 

*A2, •/' ,A4, '   LAT.  ' ,12, 

*«-■ ,F5.2,A1/26X, 'DEPTH 


«,A6,'  DATE 


•-« ,F5.2,A1, ' 


LONG. 


-'•    f 


MM. 
NB  =  2 
KK  =  0 
K  =  0 

IF(FRWTR.NE 
PHIA=-12.0 
GO  TO  1 


LENGTH  CF  CORE  ' ,F7.0, 


FROM  TOP  CF  CORE  • , F7 


MM, 


) 


,A3, 
'  ,A2 
'  ,13, 
0, 


/ 


0.  )  N2  =  3 


SAVE  PHI  6  FRACTION  WT .  FROM  EACH  DETAIL  CARD 
50  KK=KK+1 


56 


K  =  K+1 

PHI(K) =  PHIR 
FRWT(K)=FRWTP 
IF(NE.GE.8)  GO  TC  100 

SET  FLAG  (N2=2)  IF  NEW  PHI  LESS  THAN  LAST  PHI 

IF(PHIA.GE.PHI R)  N2=2 
PHIA=PHIR 
V  =  M  +  1 

IF(NINOT-l)  1,1,2 
ICO  M  =  K 

DETERMINE  FRACTION  PERCENTS  AND  ACCUMULATED  % 

DO  67  K=1,M 

PRCT(K) =(FRWT(K)/PAWT)*100. 

SUMPC=SUMPC+PRCT(K) 

£CPC(K)=SUMPC 
67  CCNTINUE 

CALL  TVAL(T,ACPC,M) 

IF(SUMPC.LT.99 .94.0R . SUMPC . GT . 100 .06 )  GO  TO  1500 

IF(M2-2)  103,1666,1990 
103  WRITE(6,802)  ( PHI ( J )  ,FRWT ( J) , PRCT ( J )  , ACPC { J)  ,  J=l , KK) 
802  FCRMATt'O' ,40X, 'PHI     SAMPLE   FRACTION   ACCUMULATED', 
*/41X,'SIZE    WEIGHT   PERCENT      PERCENT «/ (40X  ,  F6  .2 , 
*F8.4,F8.2,F12.  2)  ) 

KJK=KK 

WRITE     (6,79)     PAWT 
79    FORMATt     /40X, « POST-ANALYTICAL    WEIGHT    IS',F9.4) 

SUMNL=ACPC(KK) 

NSSC=2 

NTRSK=1 

NINM=1 

NFAW=1 

JJ=0 

DC    77    1=1, KK 

IF(FRWT( I) .EQ.O. )     GC    TO    77 

JJ=JJ+1 

T{JJ)=TU) 

ACPC(JJ)=ACPC( I ) 

PHI(JJ)=PHI(I) 
77    CCNTINUE 

KK  =  JJ 

K JK=JJ 

IF(SUMNL.GE.72.0)     GO    TO    105 

NSSC=1 

WRITE (6, 809)    PHI(KK) 
839    FOPMATC  0'  ,29X,' DID    NOT     INTERPOLATE    ANY    PHI    SIZES', 
*'     BECAUSE1 /29X, 'ACCUMULATED    PERCENT    AT    »,F5.2, 
*•     DID    NOT    EXCEED'/29X,  '72    PERCENT') 

GO    TO    160 
105     IFt4.LT.KK)    GO    TO    111 

WRITE(6,815)    KK 
815    FORMATt '     «,28X,«     CNLY     ',13,'     DETAIL    CARDS     SO    ONLY', 
*'     SAND. «/30X,  '  SILT,    CLAY    RELATIONSHIPS    CALCULATED.') 

GC    TO    901 

COMPUTE    T-VALUES 

111    CALL    INTPPtACPCPHI  ,T,KK) 
160    CALL    SNSTCLtPHI , ACPC, SUMPC, KJK) 
CALL    CTIFW 

INCREMENT    COUNTERS    AND    GO    TO    NEXT    SAMPLE 

901     KfK=KKK+KK+l 
KSM=KSM+1 
1=1  +  1 
NIN0T=2 

IFtNE.NE.9)    GO    TO    9C2 


5"7 


PRINT    RESULTS    AND    MESSAGES 


950 
860 


863 


864 


861 


959 
865 


960 
1666 

8666 


1990 
8990 


15C0 

830 


1501 
831 


833 

S323 

] 

999 


WRITE(6,86 
FORMAT (« 1' 
'     FROM     ',1 

•  CARDS.'  ) 
IF(M3T.GT. 
EXMBtL ) »L= 
FCRMAT( «0' 

•  SAMPLES' 
(33X,A7,A2 
IFCMC.GT.O 
EXMC(L) ,L= 
FORMAT  ('  0' 
«  FOLLOWIN 
•DEPTH«//( 
PTT=MBT+MC 
IFIMTT.GT. 
WRITE(6,86 
FCRMAT(  '  0" 
'  FOUND  IN 
C-C  TO  960 
WRITE  (6,8 
FORMAT ( «0" 
«     ERRORS    0 

•  MORE    CAR 
RETURN 
WRITE(6,86 
FORMAT! »0« 
•BELOW  ) 
NBT=MBT+1 
CRMB(MBT)= 
ICR(MBT  )=N 
STMB(MBT)= 
EXMB(MBT)= 
GO    TO    1501 
WRITE(6,89 
FORMAT ( «0« 

•  VALUES  B 
MC=MC+1 
CRMC(MC)=C 
JCR(MC) =NC 
STMC(MC)=S 
EXMC(MC)=E 
GC  TO  1501 
WRITE(6, 83 
FORMAT ( « 0« 
'EQUAL  POS 
F8.3,  '  .  C 
WRITE  (6,83 
FORMAT! ' 0' 
/41X,«  SIZ 
WRITE (6,83 
J=1,KK) 
FORMAT {' 
WRITE(6, 
FORMATC  ' 
) 

MED=MED+1 
GO  TO  901 
RETURN 

END 


0)    KSM,KKK 

,30X,' THIS    BATCH    0 

4,'     SAMPLES'/30X,« 


F    CARDS    CONTAINED    DATA", 
FOR    A    TOTAL    OF    ' ,14, 


0)  WR 
1,MBT 
,29X, 
/31X, 
t  A3,  A 
)  WRI 
1 ,  MC  ) 
,29X, 
G  ST  A 
23X,A 
+  MED 
0)  GO 
1) 

,30X, 
THI  S 


I TE (6,863)     (CRMB(L)  ,ICR(L) ,STM8(L), 

) 

•CARDS    OUT    OF 

•CRUISE' ,5X, ' 

2)  ) 

TE(6,864)     (CRMC(L)     , JCR( L ) , STMC ( L ) , 


ORDER    ON    THE    FOLLOWING', 
SAMPLE     • ,5X, • DEPTH'//, 


•NO  ZERO  PERCENT  CARDS  ON  THE', 
T  IONS' /31X, ' CRUISE,  5X , » SAMPLE  • 
7,A2,A3,A2) ) 


5X, 


TO  959 

• CONGRATULATI 
BATCH  OF  CAR 


65)     MTT 

,29X, • SORRY  OLD  CH 
N  THE  DATA'/30X,F0 
EFUL' ) 


ONS     NO  ERRORS  WERE', 
DS'  ) 


AP,  BUT  YOU  MADE  • ,13, 

R  THIS  RUN.  NEXT  TIME  BE  ' 


66) 

,29X, 'CARDS 


CUT  CF  ORDER.  CHECK  VALUES', 


PERCENT  CARD. '/30X, 'CHECK'  , 


CRUZ 
CR 

STAT 
EXC 

90) 

,28X,'N0    ZERO 

ELOW. ' ) 

RUZ 
R 

TAT 
XC 

0)    PAWT 

,32X,'SUM  CF  FRACTION  WEIGHTS  DID  NOT  ', 

T  ANALYTICAL  WEIGHT  '/32X, 'WHICH  WAS 

HECK  THE  VALUES  BE 

1) 

,40X,'  PHI   FRACTI 

E      WEIGHT  PERCE 

3)     (PHI ( J) , FRWT(J  ) 


LOW    FOR    ERRORS.'  ) 

ON       FRACTION       4CCUM.       T-« 
iNT  PRCT       VALUE') 

,PRCT( J) ,ACPC( J)  ,T(J  ), 


C« 
83 
0' 


,40X,F5.2,F9 
33)     SUMWT 


3,F9, 


2,F8.2,F7.3) 

'7        -  i''\i    i  K;!-:    WEIGHTS    =    ',F8.3,'     GRAMS' 


58 


c 

C      INTERPOLATION  SUBPROGRAM 
C 


1101 
1102 


151 
152 


153 
155 


157 

156 
159 
168 


169 

11S7 

154 
195 

196 

197 

199 


SUEROUT 
DIMENSI 
COMMON 
REAL 
.XTT( 3)/ 
PHIS  (6, 
PHIS  (6, 
PHIS(7, 
PHIS(7, 
PHIS(8, 
PHIS(8, 
NINT=0 
IF(NINT 
PHIS(NI 
PHIS  (NI 
NINT=NI 
XPC=PC{ 
XT=XTT( 
DO  151 
IFUCPC 
CONTINU 
YYA=PHI 
YYL=PHI 
GO  TO  2 
IF(XPC- 
IFCSUMN 
NINT=9 
GO  TO  1 
IF (XPC. 
IFISUMN 
IF(L-KK 
IF(SUMN 
IF(SUMN 
NINT=10 
GO  TO  1 
IFCSUMN 
MNT  =  11 
XT=XTT( 
GC  TO  1 
IFtL.LE 
LA=2 
GO  TO 
LA=1 
GO  TO 
LA  =  3 
LS=L-LA 
X=XT 
X1=T(LS 
X2=T(LS 
X3=T(LS 
X4=T(LS 
Y1=PHI ( 
Y2=PHI ( 
Y3=PHI ( 
Y4=PHI( 


INE  INTRP(ACPC,PHI,T,KK  ) 

ON  ACPC(  1)  ,PHI  (1)  ,T(1  ) 

/BLK1/PHIS(8,2) ,NINM,NFAW,NSSC,SUMNL,NTRSK 

PC(8)/l.i5.fl6.t25.v50.f 75.t 84.*95./t 

-2. 3 2 5, -1.645, -.995, -.674, 0.0,. 674.. 995, 1.645/ 

1  )  =  99.99 

2)=99.99 

1)=99.99 

2)=99.99 

1)=99.99 

2)=99.99 


.EQ. 
NT,  1 
NT,  2 
NT+1 
NINT 

nint: 

L  =  l,l 

.(L)-XPCJ  151,152,153 

E 

(L) 

(L) 

10 

75.0 

L.GE, 


0)  GO  TO  1102 

)=YYA 

)=YYL 

) 
} 

KK 


)  154,155,157 

,75.0)  GO  TO  156 


197 

GT.84.0)  GO  TO  159 
L.LT.84.)  GO  TO  168 
J  195,197,197 
L-95.0)  169,156,156 
L.LT.81.0)  GO  TO  1475 

197 

L.LT.92.0)  GO  TO  1476 

NINT-3) 

97 

.2)  GO  TO  196 

199 

199 


) 

+  1) 
+  2) 
+  3) 
LS) 
LS+1) 
LS+2) 
LS+3) 


AITKENS    FOUR    POINT    INTERPOLATION 

P12  =  ( (Y1*(X2-X)  )-(Y2*(Xl-X)  )  )/(X2-Xl ) 
P13=( (Yl*( X3-X) )-( Y3MX1-X) ) ) /(X3-X1 J 
P14=(  (Yl*(  X4-X)  )-(Y4*(Xl-X)))/(X4-Xl) 
P123=((P12*(X3-X) )-(P13*(X2-X J )  )/ (X3-X2) 
P124=( (P12M X4-X) )-(P14*(X2-X) ) ) /(X4-X2) 
YYA=(  (  P123*(X4-X  )  )-(  P124*(X3-X)  )  J/IX4-X3) 


59 


LINEAR    INTERPOLATION 


213 


803 


1275 


1300 
1301 


1475 


1476 


1400 
1401 


X11=T(L- 
X22=T(L) 
Y11=PHI ( 
Y22=PHI { 
YYL=(X-X 
IF(L.LE. 
IFCNINT. 
INI=NINT 
PHIS (NIN 
PHIS(NIN 
WRITE (6, 
FORMAT {• 
/32X,' 
,7X, '84' 
•     LINEAR 
RETURN 
NM3=NINT 
PHIS  (NM3 
PHISINM3 
IPH=NM3 
IPER=PC ( 
IF(NINT. 
WRITE     (6 
FORMAT! ' 
'     LEVEL 
'25' ,7X, 
30X,F5.2 
RETURN 
IPER=84 
NINM=2 
GO    TO    14 
NFAW=2 
IPEP.=  95 
GO    TO    14 
WRITE     (6 
FORMAT { • 
«     LEVEL 
'25« ,7X, 
30X,F5.2 
RETURN 
END 


1) 

L-l) 

L) 

11)*(Y22-Y 

2)     YYA=YYL 

LT.8)     GO    T 

-7 

T,1)=YYA 

T,2)=YYL 

803)     ((PHI 

0'  ,29X, «  PH 

1    (C) 

,7X, »95'/2 

•  ) 

-3 
,1)=YYA 


11)/(X22-X11)+Y11 
C    1101 


S(II,K) ,11= 
I     SIZES    AT 
5  16' 

7X,8F9.2,  • 


1,8)  ,K=1,2) 

PERCENT    LEVEL 


OF 


PT. 


(M) 


75 


,W  =  2 

:R, ((PHISCIIfK), 11=1,8       ),K=1,2) 

II     SIZES    AT    PERCENT    LEVEL    OF  (• 


,2)=YYL 

NM3) 

EQ.ll)     NFA 

,1301)     IPE 

0'  ,29X,  '  PH 

EXTRAPOLATED)  «/32X,l(C) • ,6X, • 5«  ,7X  ,  • 16 «  ,7X 

•50(M)        75 

,7F9.2,'     L 


00 


12, 


' ,7X,«84'/30X,F5.2,7F9.2, ' 
INEAR* ) 


4    PT 


T 


)tK=l,2) 


00 

,1401)     IPER,  UPHISU  I,K)  ,  11  =  1 

0,,29X,«PHI     SIZES    AT    PERCENT    LEVEL    OF     (',12, 

EXTRAPOLATED.) • /32X , 1 ( C ) • , 6X, • 5 ' , 7X , • 16 

•50(M)        75 

,7F9.2,'     L 


•  ,7X,'84,/30X,F5.2 ,7F9 
INEAR' ) 


,7X, 
PT.'/, 


C 

CO*  *t-  *X*  O-*  *A-  *JL*  *JL  *JU  O-h   OC  %1*  t£t  *i*  *J*  ■^  *V  ^*  *Jv  *•■*'  *^<  ^*  *^  ^"  *^t  *J*   *i^  ^*  *V  *J*  *&*  O*  ^*  ^^  ***  ^<  ^V  *A*  *J*  V'  ***  ^*  ^V 
T*  *T~   -  '   '"    '  *T*  -i-  -!*■    -^   ^jK  .-(»   *A   *fi  T  "   ""*  ^  "--  '^  *V*   T   "V   "V   t*    *Tt*   T'rTTTn*T'T'ilT1'''l'T   *T   *P 

C 
C 

c 
c 


-J-  «X.  -J-  Oi 


SUBPROGRAM  WHICH  I 
HAND,  THE  GRAIN  SI 
PERCENTAGE. 


*'  ■■  -J-  ~i.    »(,  »».  J,  „'-  O,  ,U  J.  .J,  J,  .1 .  ,< - 

—  -I*  -V  -V  '(*  T  1*  1>  T  1-  T  T  T  1* 


S    THE    EQUIVALENT    OF    PLOTTING,     BY 
ZE    AGAINST    THE    ACCUMULATED 


/-  -j-  -j;  o.  j.  u*  «^-  -  •■';  ^  ^-  -^  y:  »■;  J:  i  -1;  J-  i';  ;';  ^';  -'■  ^  ^  i  ^t  V:  J;  ^:  -';  "'r  :'c  V;  V'  ^J.  J;  ;';*-'-■';  5';  V-  V;  i'.   r1;  ,:?'-'  51*  -V  ;'*'  -1; 


^L*  J*  *4*  */*  J*  OL*  ^*  *JL*  *-Q  *1*  «»A»  »*0  ^U  *X*  *** 


T*  ^t~  *r"  1*  "»'  ' 


SUBROUTINE    TVA L ( T , ACPC  ,M ) 

CIMENSICN    T(  1) , ACPC( 1) 

REAL    TBLPC(87)  / 0 . 0 , 01 . 20 ,03 . 59 , 04 

1  10.  64, ll. 79, 12. S3, 14. 06, 15.  17,  16 

*  20.54, 

2  21.57,22.57,23.57,24.54,25.49,26 

*  ,29.95, 

3  30.78,31.59,32.38,33.15,33.89,34 

5};  "2  "7       29 

4  3  7.90,3  8.49,3  9.07,39.62,40.  15,40 

*  ,42.51, 

5  42.92, 42.32,43.70,44. 06,44.  41,44 
-     ,45.91, 

6  46.  16, 46. 41, 46. 64, 46. 86, 47.  06, 47 

*  ,47.93, 

7  4  8. 08,4  8.21,4  8.34,4  8.46,48.  5  7,48 

*  ,49.04, 

8  49.11,49.20,49.31,49.40,49.60,49 


.78,05.96,07.14, 08.32 
.28,17.36  ,15.44, 19.50 
.42,27.34  ,28.23,29.10 
.61 ,35.31  ,35.99,36.65 
.66, 41.15,41.62,42.07 
.74,45.05,45.35,45.64 
.26,47.44,47.61 ,47.78 
.68  ,48.78, 48. 87,48. 96 
.70, 49. 80,49.90, 50. 0/ 


60 


REAL    TB LT( 87 )/ 0.00, 0.03, 0.09, 0.12, 0.1 5, 0.18, 0.2 1,0. 24, 

*  0.27,0.30, 

1  0.33,0.36,0.39,0.42,0.45,0.48,0.51,0.5  4,0.57,0.60, 

*  0.63,0.66, 

2  0.69, 0.7 2, 0.75,0. 78, 0.81, 0.84, 0.87, 0.90, 0.93, 0.96, 

sfc       fl  o  Q  v     n  9 

3  1  *  05,'  1 '.  C8,'  1.11,  1.14, 1.17, 1.20, 1.23,  1.26, 1.29, 1.32, 

*  1.35,1.38, 

4  1.41, 1.44, 1.47, 1.50, 1.53, 1.56, 1.59,1. 62, 1.65,  1.68, 

*  1.71,1.74, 

5  1.77, 1.80, 1.83, 1.86, 1.89, 1.92, 1.95, 1.98, 2. 01,  2.  04, 

*  2.07,2.10, 

6  2. 13, 2. 16, 2. 19, 2. 22, 2. 25, 2. 28, 2. 31, 2. 34, 2. 37, 2. 41, 

*  2.46,2,51, 

7  2.65,2.75,2.88,3.08,4.09/ 

C  CALCULATE    T-VALUE 

DO    68    K=1,M 

IF  (ACPC(K)  .LT. 130.0)    GO    TO    62 

61  T(K)=4.C9 
GO    TO    68 

62  DLPC=ACPC( K)-53. 
IF(DLPC)    64,63,65 

63  T(K)=0.0 
GO    TO    68 

64  GMPC=-DLPC 
GC    TO    66 

65  GMPC=DLPC 

66  DO    69    L=l,87 
IF(TBLPC(L)-GMPC )    69,70,71 

69  CONTINUE 

70  TCALC=TBLT(L) 
GC    TO    72 

71  TCALC=(GMPC-TBLPC{L-1))*(TBLT(L)-TBLT(L-1)  )/ 

*  TBLPC(L )-TBLPC(L-l) )+TBLT(L-l) 
IF(TCALC.GT.4.09  )    GO    TO    61 

72  IF(DLPC.GE.O.O)     GO    TO    74 
TCALC=-TCALC 

74    T(K)=TCALC 
68    CONTINUE 
RETURN 
END 
C 

c 

C  SUBPROGRAM    FCR    CCMPUTING    S ANC-S ILT-CLAY    RELATIONSHIPS 

C 

C^V  -j-    -)^  «x-  »/»  *x*  -J-  «J-  .,',  *'-■  *&*  -Jr  »**  **••  **■*  ^V  *V  "V  "Ji*  "V  *&*  'J'  *ff  **•  *t*  •J*  ^M  rfc*  *A*  *'*  W*  *fr  J*  ^V  "J*  ***  *■*•  *>V  ^*  sV  *^#  *A"  **•■  ***  *}*  <4*  ■*•*  *fr  •>iV  M*  *Jr  ^*  *V  »k  *V  J*  V*  */■*  V-* 
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SANDY 
'/ 


SUBROUTINE    SNSTCL(PHI ,  ACPC , SUMPC , K JK ) 

DIMENSION    PHK1)  ,T  ( 1 )  ,  ACPC  (1 ) 

REAL*8    CLASS(9)/'     SAND  SILT  CLAY 

1,'SILTY       CLAYEY  -SILT         -CLAY 

INTEGER    SU BS (3 3) /l, 7, 7, 2, 7, 7, 3, 7, 7, 1,8, 9, 5, 1,7, 4, 2, 7, 6 
*, 1,7,6, 2,7, 
14,3,7,5,3,7,7,  7,7/ 

NI=0 

RATI0=99S9.99 

GRSN=0.0 

S4ND=0. 0 

SILT=0.0 

CLAY=0.0 

SANDP=0.0 

IF(PHI (KJK).LT.-l.O)    GO    TO    380 

CO    302    KG=1,KJK 

IF (ABS(DHI (KG) +1.0)  .LE..0001)     GO    TO    303 

IF(PHI(KG)+1.0)     202,303,307 

302  CCNTINUE 

303  GRSN=ACPC( KG) 

3C7     IF(PHI (KJK) .LT.4.0)    GO    TO    312 
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CC    309    KS= 

IF(ABS(PHI 

IF(PHKKS) 
3  09    CCNTINUE 
310    SAND=ACPC( 

SANDP=SAND 

FMUD=SUMPC 

IF(FMUD.NE 

GO  TO  313 

SAND=SUMPC-GRSN 

SANDP=SAND 


1,KJK 

(KS)-4.0) .LE..0001)  GO  TO  310 

-4.0)  309,310,313 

KS) 

-GRSN 
-SAND 
.0.)     RATIO=SAND/FMUD 


312 

313 
314 


GO    TO    380 

IF(ABS(PHI 

IF(PHI(KJK 

CO    315    KSL 

IF(ABS(PHI 

IFtPHHKSL 

315    CCNTINUE 

3  16  SILT=ACPC( 
CLAY=SUMPC 
GO    TO    320 

317    SILT=SUMPC 

320    NI=1 

IFISAND.GE 

NI  =  2 

IF(SILT.GE 

NI=3 

IFtCLAY.GE 

NI=4 

IF(SANDeGE 

*  GO    TO    3 
N4  =  l 

IF(CLAY/SI 
IF(SAND/SI 
GO    TO     (334 

334  IF(CLAY/SA 
GC    TO     (340 

336  M=5 
GO  TO 

337  NI=6 
GC  TO 

338  NI=7 
GO  TO 

339  NI=8 
GO  TO 

340  M=9 
GO  TO 

341  NI=10 

380    WPITE(6,80 

804    FORMATCO" 

1'       GRAVEL 

2  ,' SAND/MUD 

NN=3*(NI-1 

N1=SUBS(NN 

N2=SUBS(NN 

N3=SUBS(NN 

WRITE(6,87 

*  CLASS(M2) 
1    CLASS(N3) 

87    FQRMATC     • 
RETURN 
END 


(KJK)-8.0) .LE..0001)     GC    TO    314 

).LT.8.0)    GO    TO    317 

=1,KJK 

(KSD-8.0)  .LE..0001)    GO    TO    316 

).GE.8.0)    GC    TO    316 

KSU-SAND 

-SAND-SILT 

-SAND 

.75.0)     GC  TO  380 

.75.0)    GO  TO  380 

.75.0)    GO  TO  380 

.2  0. O.CR. SILT. GE. 20.0. OR. CLAY. GE. 20.0) 
80 

LT.LT.l.)     N4=2 

LT.LT.l.)    GO    TO    (341, 334), N4 

,336) ,N4 

ND.LT.l.)     GO    TO     (338, 337), N4 

,339),N4 


380 
380 
380 
380 
380 


4) 

,28X,'  SAND,  SILT,  CLAY,  RELATI CNSHI PS« /30X , 
SAND      SILT      CLAY     TOTAL    « 
TYPE' ) 

)  +  l 

) 

+  1) 

+  2) 

)    GRSN,SANDP,SILT, CL AY, SUMPC , R AT  10 , CLA SS( Nl ) 

, 

,29X,F8.2,5F9.2,3X,3A6) 
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C 

C      SUBPROGRAM  FCR  COMPUTING  TRASK,  INMAN,  AND  FOLK  AND 

C      WARD  STATISTICS. 

C 

SUBROUTINE  CTIFW 

CCM.MCN    /BLK1 /PHIS  (8,2)  , N INM, NFAW, NSSC ,SUMNL,NTRSK 

REAL*8    ZMEN(20) /' VERY    COARSE     SANDCCARSE    SAND  MECIU 

IN    SAND  FINE    SAND  VERY    FINE    SAND       COARSE     SILT 

2  MEDIUM    SILT  FINE    SILT  VERY    FINE    SILT 

3CLAY  «/ 

REAL*8    DEV  1(18 J/'VERY    WELL    SORTED  WELL    SORTED 

MODERATELY    SORTED  POORLY    SCRTED 

VERY    POORLY    SCRTED  EXTREMELY     POORLY    SORTED'/ 

REAL*8    SKEW(15 )/'VERY    NEGATIVELY    SKEWED       NEGATIVELY    SK 
1EWED  NEARLY    SYMMETRICAL  POSITIVELY    SKEWED 

2  VERY    POSITIVELY    SKEWED        •/ 

REAL*8    KURT(15 ) / • PL AYTKURT IC  MESOKUPTIC 

1  LEPTGKURTIC  VERY    LEPTCKURTIC 

2  EXTREMELY    LEPTOKURTIC  •/ 

REAL    TLFC(7)/0.0,0.35,J.50,1. ,  2  .  ,4. ,99. 99/ ,TLFS ( 6 )/-3. 

&       —        ^       —         1 

1  13  t9!99/>TL1=K(6)/0.f  .9  ,  1.  11,  1 .5,  3.  ,99.99/ 

REAL    HEAD(6)/«4     PT.     PHIS       LINEAR    PHIS    •/ 
INTEGER     ITLFSC(6)/2*1,2,3,4,5/ 
DC    1000    LL=1,2 
IF(NSSC.EQ.l)     GC    TC    400 

C  CALCULATE    TRASK    VALUES 

PFI5=PHIS(2,LL ) 

PHI16=PHIS  (3,LL) 

PHI25=PHIS(4,LL) 

PHI50=PHIS( 5,LL) 

PFI75=PHIS(6,LL) 

PHI84=PHIS(7,LL) 

PHI95=PHIS(8,LL) 

Gl=2.** (-PHI25) 

Q2=2.**(-PHI50) 

Q3=2.**(-PHI75 ) 

SC=SQRT(Q1/G3) 

FLGS0=AL0G10(S0) 

SKG=SQRT((Q1*Q3)/(Q2*Q2) ) 

LLL=3*(LL-1)+1 

LLT=LLL+2 

WRITE (6, 80 5}  (HEAD(LLK) ,  L LK  =  L LL , LLT ) , Q 1 , Q2 , Q3 , SO,  FLGSO, 
*2SX, «    Ql        Q2        Q3        SC      LCG  SO 
*SKG,/30X,F6.3,4F9.3,F8.2) 

IF(NTRSK.GT.l)  GC  TC  1000 

C      CALCULATE  INMAN  VALUES 

IF(NINM.GT.l)  GO  TO  505 
FIVD=( PHI16+PHI84)/2.0 
FIDV=(PHI84-PHI16)/2.0 

FISK=(FIMD-PHI 50) /FIDV 

IF(PHI95.EQ.  99.99)  GO  TO  504 

F2SK=( (  (PHI95  +  PHI5 )/2. ) -PH I  5 0 ) / F I DV 

FIKU=( ( (PHI95-PHI5J/2. )-FIDV  )/FIDV 

W«ITE(6,806  )  PHI50,eiMD,FIDV,FISK,F2SK,FIKU 
806  FORMAT(  '  0'  ,29X,« INVAN  VALUE S ' /32X , ' ME C I  AN     MEAN',5X, 
l'CEV.      SKEWED.     2ND  SKEW.      KURT.'/32X, 
2F5.2,3F9.2,2F11.2) 

GC  TO  600 

504  WRITE(6,92)  PH I  50, F IMD , F I D V , F I SK 

l'ND  SKEWNESS  AND  KURTOS IS )• /30X ,« MEDIAN     MEAN       ' 
1,'DEV.      SKE W. «/30X,F5.2 ,3F9.2) 
GC  TO  600 

505  WRITE(6,<53) 
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93  FORMAK  •  0' ,28X,  •  INMAN  PLUS  FOLK  AND  WARD  VALUES  NOT  « 
l'CALCULATED  BEC AUS E • /29X ,  «  NEXT  TO  LAST  ACCUMULATED  « 
2PERCENT  WAS  LESS  THAN  84') 

GO  TO  1000 
400  WRITE{6,91) 
91  FORMAK  '0'  ,28X,«  NCT  ABLE  TO  CALCULATE  TPASK,  INMAN,  « 
*,'  OR  FCLK  AND  WARC/29X, 
*'  VALUES  BECAUSE  NEXT  TO  LAST  ACCUMULATED  PERCENT  DID' 
*,'  N0T'/29X,«  EXCEED  72') 
GO  TO  1000 

CALCULATE  FOLK  AND  WARD  VALUES 

600  IF(NFAW.EQ.l)  GO  TC  601 
WRITE  (6,94) 

94  FORMAT( '0' ,29X, 'COULD  NOT  CALCULATE  FOLK  AND  WARD' 

*■  VALUES  BECAUSE  NEXT' /30X , • TO  LAST  ACCUMULATED  PERCEN 
*,«NT  DID  NOT  EXCEED  92.') 

CCMPUTE  MEAN  AND  CETERMINE  CATEGORY 

6  01  FMZ=(PHI16+PHI50+PHI84)/3.0 
DC  300  1=1,10 
ZVAL=I-1 
IFIFMZ.LT.ZVAL )  GO  TO  301 

300  CONTINUE 
1=10 

301  Il=2*(  I-D+l 
12=11+1 

WRITE     (6,8C7)     FMZ , ( ZMEN ( 1 1 )  , 1 1 =1 1 , I  2 ) 
8C7    FORMATCO' ,29X, 'FOLK    AND    WARD    VALUE  « /31X  ,  •  MEAN  •  ,  F  12.  2  , 
14X,2A8) 

CGMPUTE    DEVIATION    AND    DETERMINE    CATEGORY 

FDEV=(PHI84-PHI16) /4.+  ( PH I 95-PHI 5 ) /6.  6 

DO    604    L=l,7 

IF(TLFD( L) .GE.FDEV)     GO    TO    605 

604  CONTINUE 

605  IFDTL=L-1 

IF( IFDTL.EQ.O)  IFDTL=1 
Il=3*(  IFDTL-D+1 
12=11+2 

WRITE  (6,131)  FDEV,  (CEVI  (II)  ,11  =  11,12) 
131  FORM AT (3 IX, 'DEVIATION' ,  F7. 2, 4X , 3A 8 ) 

COMPUTE  SKEWNESS  AND  DETERMINE  CATEGORY 

FSK=(PHU6  +  PHI84-2.0*PHI50)/( 2. *( PHI 84-PH I  16)  )  + 
1  ((PHI  5  +  PHI95)-(2.*PHI50)  )/ ( 2.* ( PHI95-PH I  5 ) ) 

DO  6  08  L=l,6 

IF(TLFS (L) .GE.FSK)  GO  TO  609 
6C8  CONTINUE 
6C9  IFSKTL=ITLFSC(L) 

I1=3=M  IFSKTL-D+1 

12=11+2 

WRITE{6,142)  FSK,  (  SKEW(  II)  ,  11  =  11,  12) 
142  FCRMAT(31X,' SK EWNES S ' , F8 .2 ,4X , 3 A8 ) 

;      COMPUTE  KURTCSIS  AND  DETERMINE  CATEGORY 

FKG=(  PHI  95- PHI  5)  /  (  2.44*  ( PHI75-iPH  125 )  ) 

DO  612  L=l,6 

IF(TLFK(L) .GE. FKG)  GO  TO  613 

612  CONTINUE 

613  IFKTL=L-1 
IF(IFKTL.EG.O)     IFKTL=1 
Il=3*(  IFKTL-D+1 
12=11+2 

WRITE (6, 16 2)     FKG, (KURT (II) ,11=11,12) 
162    FORMAK 3 IX,  'KURTOSI S'  ,  F8  .  2  ,  4X  ,  3A 8 ) 
1000    CONTINUE 


6U 


RETURN 
END 
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